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ABSTRACT 

Context. The self-enrichment scenario for globular clusters (GC) requires large amounts of residual gas after the initial formation 
of the first stellar generation. Recently, Krause et al. (2012) found that supemovae may not be able to expel that gas, as required 
to explain their present day gas-free state, and suggested that a sudden accretion on to the dark remnants, at a stage when type II 
supemovae have ceased, may plausibly lead to fast gas expulsion. 

Aims. Here, we explore the consequences of these results for the self-enrichment scenario via fast rotating massive stars (FRMS). 
Methods. We analyse the interaction of FRMS with the intra-cluster medium (ICM), in particular where, when and how the second 
generation of stars may form. From the results, we develop a timeline of the first x 40 Myr of GC evolution. 

Results. The results of Paper I imply three phases during which the ICM is in a fundamentally different state, namely the wind 
bubble phase (lasting 3.5 to 8.8 Myr), the supernova phase (lasting 26.2 to 31.5 Myr), and the dark remnant accretion phase (lasting 
0.1 to 4 Myr): (i) Quickly after the first generation massive stars have formed, stellar wind bubbles compress the ICM into thin 
filaments. No stars may form in the normal way during this phase, due to the high Lyman-Werner flux density. If the first generation 
massive stars have however equatorial ejections, as we proposed in the FRMS scenario, accretion may resume in the shadow of the 
equatorial ejecta. The second generation stars may then form due to gravitational instability in these discs that are fed by both the 
FRMS ejecta and pristine gas. (ii) In the supernova phase the ICM develops strong turbulence, with characteristic velocities below 
the escape velocity. The gas does not accrete neither on to the stars nor on to the dark remnants in this phase due to the high gas 
velocities. The strong mass loss associated with the transformation of the FRMS into dark remnants then leads to the removal of the 
second generation stars from the immediate vicinity of the dark remnants, (iii) When the supemovae have ceased, turbulence decays 
quickly, and the gas can once more accrete, now on to the dark remnants. As discussed in Paper I this may release sufficient energy to 
unbind the gas, and may happen fast enough so that a large fraction of less tightly bound first generation stars are lost. 
Conclusions. The detailed study of the FRMS scenario for the self-enrichment of GCs reveals the important role of the physics of 
the ICM in our understanding of the formation and early evolution of GCs. Depending on the level of mass segregation, this sets 
constraints on the orbital properties of the stars, in particular high orbital eccentricities, which should have implications on the GC 
formation scenario. 

Key words. (Galaxy:) globular clusters: general - ISM: bubbles - ISM: jets and outflows 

1 . Introduction (SNe) ejecta have to be removed in order to keep the constancy 

of [C+N+O] and the mono-metallicity observed in most GCs. 

H ■ Recent detailed spectroscopic and deep photometiic studies 
. . ." have lead to the conclusion that Galactic globular clusters (GCs) Two main models have been developed in the literature, 
host multiple generations of stars; those differ in their chemi- according to the nature of the first generation stars poten- 
cal abundance properties and in some cases in their membership tially responsible for the GC self-enrichment (hereafter the stel- 
to multimoda l sequ ences in the colour-i nagnitude diagrani (e.g. lar polluters), namely the fast-rotating-massive-stars (FRMS) 
iGratton et al.ll2004l 120121: iPiottol 120091; ICharbonnellbOlOt for or the asymptotic-giant-branch stars (AGB). Both the FRMS 
recent reviews). The self-enrichment scenario to explain the and the AGB scenarios face a problem in the mass budget 
observations is Ukely complex. It requires that the products of between the amount of matter provided by the polluter stars 
hot hydrogen-burning ejected by fast evolving first generation and the mass locked today in first and second generation low- 
stars, remain within the GC and are recycled in second gener- mass stars (about 30 vs 70 % respectively in most GCs, and 
ation stars to explain e.g. the ubiquitous O-Na anti-correlation; up to 50-50 in few cases; for more details see ICarretta et al.l 
but on the other hand, helium-burning products and supernova l2009l) . Solutions require either a top-heavy IMF for the first stel- 

lar generation do^ntona & CelIo] 2004; Bekk i & Morris 200^1 

iPrantzos & Charbonnell 120061 : iDecressin et al.l l2007al) . or sub- 

* E-mail: Martin.Krause@unige.ch stantial loss of first generation low-mass stars from the GCs 
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which must then have been much more massive initially than 
today ('Decressin et al."2007a'; 'Bekki et al."2007';'D'Ercole et alj 
[2008; Decressin et al. 2010; Schaerer & Charbonnel 2011 ), the 
second option being currently favoured. In the FRMS sce- 
nario, the hydrogen-burning ashes would be ejected from their 
fast rotating parent stars due to repeated super-critical rota- 
tion. These ejecta are thought to mix with pristine gas and then 
form the second generation stars in the immediate surround- 
ings of the massive star polluters. (iPrantzos & Charbonnell2006t 
iDecressin et al.ll2b07allbh . Initial or very early mass segregation 
together with quick gas expulsion after the formation of the sec- 
ond stellar population, are invoked to loose preferentially first 
generation low-mass stars and ret ain the observed propo rtion of 
first and second generation stars (IDecressin et al.luOlOl ). In the 
asymptotic-giant-branch (AGB) scenario, the pristine gas left af- 
ter the formation of the first stellar generation is supposed first to 
be expelled together with the unwanted SNe ejecta and the bulk 
of first generation low-mass stars. Later when the SNe would 
have ceased, the slow winds of AGB stars would accrete in a 
cooling flow towards the GC centre, where they would form the 
second generation of stars after dilution with re-collected pris- 



tine gas whose sources are still uncertain ( 


D'Ercole et al.ll2008l 


l20iai201lLl2012l;IConrov & SDer2ei2011 


3. 



In lKrause et alj ( 1201 2l Paper I) we showed that gas expulsion 
via SNe, which has long been the prevailing paradigm to change 
the GC potential well and induce the loss of first generation low- 
mass stars, does not work in GCs. The reason is that, while the 
energy produced by SNe usually exceeds the binding energy, it is 
not delivered fast enough to avoid the Rayleigh-Taylor instability 
of the escaping shell. We have also shown that a more powerful 
event, such as the energy released by accretion on to the dark 
remnants, may lead to successful gas expulsion. These results 
suggest three distinct chronological phases in typical GCs: Wind 
bubbles are created during the whole period from about to 8.8 
Myr on our global time scale. Core collapse to black holes and 
neutron stars occurs from about 3.5 Myr until 35 Myr. Where 
black holes are formed, the collapse may not be accompanied by 
a strong or even any momentum ejection. If we therefore con- 
sider that only stars with initial masses below 25 M© will give 
birth to an energetic supernova event then energy ejection by su- 
pernovae will occur only between 8.8 and 35 Myrs on the global 
timescale. Finally, the dark remnants are activated and expel the 
gas as well as the majority of the first generation stars. 

Here, we explore the implications of these findings for the 
formation of the second generation stars in the context of the 
FRMS scenario and develop a detailed timeline for the first 
R! 40 Myrs in the lifetime of a typical GC. We describe the ba- 
sic model setup in Sect.|2] The wind bubble, supernova and dark 
remnant accretion phases are desciibed in Sects. |3] Hand|5] re- 
spectively. We discuss our findings in Sect.|6]and summarise and 
conclude in Sect. [T] In order to describe the sequence of events, 
we refer to a global timescale throughout this article. The global 
clock is set to zero at the coeval birth of the first generation of 
stars. 



Dilution of the AGB ejecta with re-collected gas with pristine 
chemical composition is mand atory since all current AGB models (e.g. 
[Karakas & Lattanzio 2007; V entura & D'A ntona 2009; Siess 20l3) pre- 
dict that O and Na yields should actually correlate, at odds with the 
observed abundance anticorrelations. Although several scenarios have 
been proposed, the gathering mecha nism of the required pristine matter 
is still uncertain (for discussion see iD'Ercole et al.l20To[|201 ih . 



Table 1. Parameters of model cluster 



Parameter 


Value 


Description 


Mot 


9x IO^'Mq 


Total mass 


n/2 


3 pc 


Half-mass radius 




1/3 


Star formation efficiency 


PO 


10'';7!p cm"-' 


Average gas density 



2. Basic model setup 

The model presented here gen erally follows the ideas outlined in 
IDecressin et al ] (l2007bll201(ih . We consider an initial dense gas 
cloud of Mtot = 9 X IO^Mq with a half-mass radius of ri/2 = 3 pc 
which forms the first generation stars at an efficiency of = 1/3 
and according to a Salpeter initial mass function (IMF) for first 
generation stars more massive than 0.8 M© and to a log-normal 
distribution for less massive long-lived stars of first and second 
generations. The parameters summarised in TablelTjare inferred 
from N-body simulations which assume a Plummer distribution 
for the spatial distribution of gas and stars (i.e., the star f ormation 
efficiency is similar at all radii; IDecressin et al ] l2010h . For our 
model cluster, we find about 5700 first generation massive stars 
between 25 and 120 Mq. Mass segregation leads to all massive 
stars to be located within ri/2- We use the stellar evolutionary 
models of FRMS at subsolar meta l licity (Z = 0.0005, [Fe/H]~ 
-1.5) presented bv lDecressin et al 1 (l2007bl) . From these models, 
we have extracted the relevant wind parameters, summarised in 
Table ID IDecressin et al. l (l2OT0l) use these parameters to describe 
the GC NGC 6752 with a present day mass of 3 x 10^ Mq, as- 
suming complete recycling of the slow wind released by FRMS 
and dilution with pristine gas to reproduce the observed Li-Na 
anti correlation. For this setup, we will now explore the different 
phases in detail. A sketch of the model is provided in Fig.[T| 

3. Hot wind bubbles 

We will now argue (Sect. 13. Il l that during the first few Myrs of 
the GC, before the first type II SNe, the fast radiative winds of 
the massive stars will create large bubbles that will probably also 
overlap, but they will probably not unite into a single superbub- 
ble and will not lift any important amount of gas out of the GC. 
This situation is depicted in Fig. [T] (second row from top). We 
also show that during this wind bubble phase (and also the super- 
nova phase), no stars should form in the normal way (Sect. 13.2b . 
the equatorial ejections of the FRMS establish a decretion disc 
close to the FRMS (Sect. 13.31 ) and that accretion may be re- 
established in the shadow of these decretion discs (Sect. 13.4b . 
The inner decretion and the outer accretion discs merge on the 
viscous timescale (Sect. 13.5b . and finally the second generation 
stars form in the merged discs (Sect. 13.6b . 

3. 1 . Bubble expansion 

To see that the wind bubbles may not unite into a single su- 
perbubble and that they will also not lift much gas out of the 
GC s potential well, we consider the energy budget in the clus- 
ter. iFrever et"aD (l2003l) derive that the mechanical power out- 
put of a star dominates the effect on the intra-cluster medium 
(ICM) rather than the effect of ionisation. The power output of 
the fast winds of a FRMS with a given mass M between 40Mo 
and 1 20Mo may be approximated by 

go(M) = 4xlO^°(M/Mo)2-^ergs"^ (1) 



2 



Martin Krause et al.: Superbubbles in globular cluster infancy II 



Total bubble volume within half-mass radius 




0.10 



Time / Myr 



Fig. 2. Combined volume of the wind bubbles of all the 5700 
massive stars (M > ISMq) in our model cluster as fraction of the 
half mass radius. The solid (dotted, dashed) line is for an energy 
injection efficiency of 77 = 0.1 (1, 0.01). 

This power law follows from interpolation of the values in 
Table [S] With a Salpeter IMF, the total wind output, which 
is dominated by the most massive stars, integrates to: 2o,tot - 
2 X 10-^*^ erg s-i = 6 X 10^^ gj-g Myr'^ 

How much of this energy is actually used to e xtend the bub- 
bles and not lost to radiation? lFrever et al.l (l2003l) derive that for 
a background density of 20 cm~^, much lower than the density 
expected in a young GC, about 90 per cent of this power is ra- 
diated away, and thus the energy efficiency is only 10%. If the 
bubble would remain self-similar, a greater energy e fficiency is 
expected, around 7 0% (Weaver etal. 19771: iFtever e t al. 2003]). 
iKrause et alJ (1201 3h study the energy efficiency of merging wind 
bubbles, and find that for small separations, there is no addi- 
tional effect from bubble merging and the efficiency factor is not 
changed. We expect that in the dense environment of a forming 
GC radiation losses may increase and the efficiency of energy 
transfer may be reduced below the value of iFrever et al.l(l2003h . 
On the other hand, the optical depth is higher, and more of the ra- 
diation energy may be captured by the gas. In the following, we 
adopt 7; = 0.1 as a working hypothesis for the energy efficiency, 
but the conclusions would be unchanged if 77 = 0.01. 

The gravitational binding energy of the gas is iigrav - 0.4(1 - 
esf)GM^/R ^ 6 X 10^^ erg dBaum^ardt et al.ll2008h . The com- 
bined wind power of all the massive stars amounts to //Qo.tot - 
6 X 10^ '(77/0.1) erg Myr~'. Thus, it is not possible for the stellar 
winds to lift any noteworthy amount of gas out of the GC on a 
relevant timescale. 

On the other hand, we may integrate the volume fraction 
occupied by the wind bubbles by integration of the standard 
wind bubble volume 47Trl ... /3 over the IMF. The bubble ra- 

bubble' 

dius rbubbie at time t is given by: 



(2) 



'"bubble = a t ' , 

\ PO / 

where a * 0.8 for self-similar expansion (IWeaver et al.lll977l) . 
and t he not necess arily self-similar thin shell approximation 
alike (lKrausell2003h . The average gas density in our model clus- 
ter IS po - cm . We include the efficiency factor 77 in 
Eq. |2] because the expansion law in this f orm is consistent with 
the si zes of observed wind bubbles (e.g. lOev & Garcfa-Segural 
|2004 and references therein). We show the volume fraction oc- 
cupied by hot bubbles in Fig.|2] The hot wind bubbles fill the GC 
very quickly on a timescale of 0. 1 Myr. 

Thus, while the wind bubbles are not able to lift the gas out of 
the potential well of the GC, they are nevertheless able to quickly 



fill the entire volume within the half-mass radius, no matter what 
we assume in detail for the energy efficiency. It follows that 
the bulk of the cluster gas is compressed into thin filaments or 
sheets, while most of the volume is occupied by hot wind bub- 
bles. Bubbles of stars of different masses in general do not have 
the same pressure, and thus many individual interfaces may be 
expected to be swept away and pushed against some other bubble 
walls (van Marie et al. 2012; Krause et al. 2013). So, overall it is 
clear that some individual bubbles will unite to form larger bub- 
bles, but there will not be one overall superbubble which would 
imply significant lift-up of gas out of the potential well of the 
GC. 



3.2. Lyman-Werner flux 



IConrov & Spergei (1201 ll) make the point that during the main 
sequence phase of the massive stars, the Lyman- Werner flux 
(912A - llOOA) is so high that molecular hydrogen is disso- 
ciated throughout the cluster, and no stars may form in the usual 
way during this time. We repeat their ana lysis here, but wit h the 
numbers from the FRMS models of iPecressin et al.l(l2007bl) . We 
have calculated the Lyman-Werner flux from these models and 
included it in table|2] A reasonable fit to the photon fluxes, which 
are roughly constant during the main sequence, as a function of 
mass of the parent star is: 



2lw = 7 X lO^^(M/MQf '^ s"' 



(3) 



IConrov & Spergei (1201 ll) estimate the radius of the photo- 
dissociation region to 



/?LW = 1-1 pc 



LW 



,1/3 



\ 1049 I 



lO^Ti 



-3\ 



1/3 



0.028Zf 



, (4) 



adapted to the metallicity of our model stars. Thus, a 120Mo star 
is already able to photo-dissociate a large fraction of our model 
GC. Integration over a Salpeter IMF gives the factor /lw, which 
describes how many times the GC volume could be photodisso- 
ciated by the ultraviolet radiation: 



/lw = 200 



noVocm" 



3\ 



P 



0.028Zp 



(5) 



So, even if the metallicity or the gas density would be higher by 
one or even two orders of magnitude, the gas would still be fully 
photodissociated. 

In addition to the effects of dust considered in the above ana- 
lysis, molecular hydrogen may also form via electrons, which 
first combine with neutral hydrogen to H~. H2 is then formed 
via the reaction: 



H-HH-^H2-He. 



(6) 



Regions where bot h, electrons and neu tral hydrogen atoms, are 
abundant are rare. iRicotti et al.l (1200 lb show that the the inter- 
face between an H II region and a neutral region is such a 
promising place, and that molecular hydrogen may form there 
in a narrow layer. [Ricotti et al. (2001) find a strong depen- 
dence on the input spectrum and a maximum column den- 
sity of log(A/H, /cm~^) = 14 - 15 of molecular hydrogen formed 
in this layer. The corresponding gas mass itself would be too 
small to form any sign ificant amount of stars in the GC con- 
text. However, iRicotti et al. (200JJ also show that the column 
is sufficient for some Lyman-Werner bands to become optically 
thick. Hence, such layers might in principal shield pockets of 
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Table 2. Properties of low metallicity fast rotating massive stars. 











T'eq.stait / 




A/onV 


A?off'7 


Qo'/ 




Mo 


Re 


Mq Myr"' 


Mq Myr"' 


Myr 


Myr 


Myr 


Myr 


lO-'^erg s"' 


10^^-1 


120 


12 


1 


66 


0.7 


3.15 


0.020 


0.080 


28 


9 


60 


10 


0.2 


14 


0.9 


4.4 


0.028 


0.067 


3.9 


1.26 


40 


8 


0.06 


14 


1.1 


5.7 


0.065 


0.250 


1.2 


0.192 



Notes. Mass of the star Photospheric radius. Mass loss rate of the steady fast wind. Mass loss rate of the slow equatorial wind. 

Global time at which the intermittent equatorial ejections start. ^-^^ Global time at which the intermittent equatorial ejections end. Duration of 
the individual slow wind episodes. Time span in between any two slow wind episodes. Total power of the steady, radiatively driven fast wind 
on the main sequence. Flux of Lyman- Werner photons (912A - 1 lOOA) at 1 Myr after reaching the main sequence; roughly constant during the 
main sequence. 



dense gas from photo-dissociation. Quantitatively, 'Ric otti et al.l 
(1200 ih find that the radius of photo-dissociation regions may 
be reduced by a factor of 1.5, corresponding to a factor of 3.4 
in volume. This would however be too small to be significant, 
since we have shown above that there are enough photons to 
ionise the entire volume about 200 times. This crude order of 
magnitude estimate stand s up to recent three-d i mensi onal radia- 
tive transfer simulations: IWolcott-Green et aP (1201 ll) show that 
passage through molecular hydrogen with a column density of 
= 5.x: X 10 '"^cm"^ attenuates the Lyman- Werner flux (within 
an accuracy of 15 per cent) by factor 



/shield (-^Hz^r) 



0.965 



■ + ■ 



0.035 



xexp[-8.5xl0Al+-«)*'-^] , 



(7) 



where ^75 = fe/(kms~') is the Doppler broadening parameter. 
If thermal broadening dominates, should be around or be- 
low unity. If their is substantial velocity shear in the filaments, 
it might be higher. A reduction of the Lyman-Wemer flux by 
a factor of 200 (in order to get the photo-dissociated volume 
fraction below unity) is reached at Iog(A^H2/cm~^) = 17 (16,18) 
for Iog(/75/kms~') - (-1,1). This confirms that the aforemen- 
tioned layers at the interface of ionised and neutral regions with 
Iog(A^HT /cm~^) - 14-15 would not provide a significant optical 
depth which could significantly alter the conclusion above. A 
combination of this effect together with strong dust absorption 
might perhaps allow for some molecular hydrogen in the dens- 
est parts of the ICM. However, the question would remain, why 
these dense pockets would not have f ormed stars right away. 



We therefore confir m the result of Conroy & Spergell(l201 It 
with the star models of iDecressin et al. I (I2007bh that molecular 
hydrogen is no longer present in significant amounts in the GC 
as soon as the massive stars reach the main sequence. From the 
similarity of the numbers for our star models at a metallicity of 
Z - 0.03Zo, with their case for Z = O.IZ©, it should follow that 
this result should not depend very much on the exact metallicity 
of the GC. Thus, at least low-mass stars do not form in the "clas- 
sical" way during this time, because the gas is efficiently kept 
at » 100 K. We note that this remains true throughout both, the 
wind bubble and the supernova phase. 

3.3. Equatorial mass ejections 

Angular momentum transfer within the FRMS causes episodic 
equatorial mass ejections of gas enriched by hydrogen-burning 
produ cts when the star reaches break-up, formin g decretion 
discs (IDecressin et all l2007bl: iKrticka et al] l20IIh around the 
stars. This process is analogous to the equatorial mass loss in 



Be-type stars (e. g. iRivinius et al. 1 1200 II: iTownsend et al.ll200l 
iHaubois et"aDl20I2h . In FRMS, the total amount of mass lost 
by this process is rather high, and it is hardly plausible that the 
ejected material may immediately leave the vicinity of the star: 
In order to lift the material out of the gravitational potential of 
the parent star, an energy supply at a rate of GMM jilRi^) is re- 
quired, if the decretion disk is first established at a characteristic 
radius Rm after being ejected. The factor of two in the denom- 
inator takes into account the rotational energy of the gas. This 
material leaves the star because of supercritical rotation. So, it 
is reasonable to assume that it will be ejected with the neces- 
sary angular momentum to be stable just outside the star. As an 
estimate, we take /?dd to be three photospheric radii (compare 
Table|2]i. For the 40 (60,120) star, we find a required energy 
injection rate of 1.4 (1.7,13)xlO^^ erg/s. The power in the fast, 
radiative wind is a few, up to a factor of ten times lesfl The 
equato rial ejections shoul d cover only a small part of the solid 
angle (iKrticka et al.ll201 Th . whereas at the same time the radia- 
tion of fast rotating stars is concentrated towards the poles (e.g. 
Maeder & Mevnet 2012). In the light of these considerations, it 
seems highly implausible that a significant part of the disc may 
leave the immediate vicinity of the parent star, though, there may 
likely be some moderate mass loss via photo-evaporation, and 
the disc may spread out somewhat via angular momentum trans- 
fer (^ticka et al. 201 L) . 

The mass lost by the equatorial mechanism averaged over 
the total ejection time, which includes the main sequence and the 
luminous blue variable stages, may be well fitted by the formula: 



eq 



M 

60Mf. 



-4.8 



Mo/ Myr, 



(8) 



for stars with mass M > 32Mq. The latter limit is simply due to 
the zero-point in Eq. (|8]l. This is however quite realistic, as the 
40 Mq model still has a high equatorial mass loss rate, whereas 
in the 25 Mq model it is close to zero. 

3.4. Accretion on to tlie discs around tlie massive stars 

3.4.1 . The discs persist around tlie FRIVIS 

The off-time between any two equatorial mass ejections is a 
few times the on-time (Table |2]i. Also the radiative energy ac- 
cumulated in the off-time, when no new mass is ejected equa- 
torially, but the radiation pressure keeps pushing, is likely not 
sufficient to unbind the entire disc from the parent star. Another 



^ Our wind power Qq refers to the terminal wind power of the fast, 
radiatively driven wind, but the additional power required to lift the 
material ejected in this wind out of the gravitational potential is less 
than 50 per cent of Qq. 
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process to consider is close encounters of other massive stars, 
which cou ld in p rinciple unbind a disc dClar ke & Pringlelll993t 
et al.lll996h . In the unfavourable prograde encounter, the 
disc would be stripped down to about half the periastron ra- 
dius. The timescale for a given massive star to have a close en- 
counter with another massive star is given by (ncrv)"', where 
n - 2Q,QQQ I {An r^^ 2>) - 192 pc"-' is the density of massive stars, 

o" = 7r(2/?dd)^ is the characteristic cross section of the accretion 
disc, and v ^ 100 km/s is the typical velocity of the massive stars. 
Using any of the values for the disc radii discussed above, the 
encounter timescale turns out to be much larger than the Hubble 
time. Thus, the discs cannot be tidally stripped. We note that, 
due to the small cross section of the FRMS discs, ram-pressure 
stripping is equally ineffective. 



Characteristic velocities for the model cluster 



200 



150 



100 



50 





escape 




circular 


- 


_ _ , central free-fall 


''^ 





3.4.2. Re-establishment of accretion 

Importantly, the shielding from the radiation of the parent star 
should allow a re-establishment of the accretion flow providing 
pristine gas (Fig. [1] second row from top, right). The Bondi 
accretion rate should set an upper limit for the mass accumula- 
tion in the disc. With the sound speed for an atomic gas con- 
sisting of 90 per cent hydrogen and 10 per cent helium by num- 
ber, cl - jk^T l(\.21m^), and a solid angle fraction /sd of the 
disc usable for ac cretion, the Bondi accretion rate evaluates to 
(iFrank et al.ll2002h : 



Mb 



ondi 



4600Mo/ Myr 
M 



X 



60Mo 



lOOK 



3/2 



(9) 



This accretion rate is most likely not sustained throughout 
the lifetime of the massive stars. One reason is that the dense 
gas should be compressed in filaments at this time (compare 
above), so that accretion may only take place when the star 
passes through a filament. We show below that, for suitable orbit 
parameters, the total path an FRMS traverses at low velocities is 
long compared to the typical bubble size. Since the dependence 
on the density is linear, we therefore do not need to take this 
into account for the present discussion, as a reduced accretion 
rate in voids would be compensated by an increased one in the 
filaments. 



3.4.3. Suppression of accretion due to stellar velocities 

If a star moves with a velocity v relative t o the local ICM, the ac - 
cretion rate is reduced by a factor ( v/cs)^ dBondi & Hovlell944t) . 
For typical virial velocities of about 100 km/s, this would lead to 
a suppression of about six orders of magnitude. Orbits of stars in 
GCs are not expected to be circular, and therefore the accretion 
rate is a strong function of the position of a given star along its 
orbit: It accretes more strongly when the velocity relative to the 
ICM is small, i.e. near the outer turning point and when at the 
same time it passes through a filament. 

We will now estimate the time averaged accretion rate for a 
given star. In a spherical poten tial, the orbits are general ly ex- 
pected to be of the rosette type dBinnev & Tremainell2008h . with 
most of the stars having only slightly negative total energies. 
Adapted to our model cluster, most of the massive stars will have 
just enough energy to reach the core radius at their outer turning 
points. As may be seen from Fig. [3] the circular velocities at the 
core radius (2.3 pc) are around 80 km s~' . If the stars would have 
no angular momentum, they would reach a similar velocity after 



r/ pc 



Fig. 3. Characteristic velocities for the model GC (compare 
Table [T]) as a function of radius. The solid line shows the lo- 
cal escape velocity. The dotted line is for the circular velocity. 
The dashed line shows the velocity a star would have obtained 
in the centre if it would start at rest at a given radius r. 



a free fall to the centre from this position. Consequently, veloci- 
ties around 1 km/s, as necessary for efficient accretion may only 
be reached for very eccentric orbits and only near their outer 
turning points. 

For this estimate, we approximate the orbits by ellipses. The 
radius r as a function of azimuth is then given by 



(10) 



1 -l-ecos(0) ' 



with the eccentricity e and the orbital size parameter p. The or- 
bital velocity follows by differentiation and may be expressed 
as 



\V(0) = vo 



1 H-2ecos0-He2 



(11) 



where = GMtot/a with the semi-major axis a. Instead of the 
eccentricity, we use in the following the ratio of the velocity at 
the outer turning point to the circular velocity at this location 
as parameter. It is given by Vaz.o/vcirc = V(l -«)/(! +e)vo/v'circ- 
The orbit-averaged Bondi accretion rate may then be found by 
integration of Eq. (|9]l over the appropriate number of orbits, n, 
multiplied by the suppression factor, (Cs/vg)^, and dividing by n 
times the orbital period T: 



Bondi 



>=M, 



Bondi 



nT Jo 



dt 



(12) 



We remind here that Mfiondi depends on density and tempera- 
ture. We have taken this out of the integral in Eq. (fT2t . be- 
cause the changes in density along the path of a given star - 
due to filaments and voids - are not correlated to the changes 
in velocity. Again, since the Bondi accretion rate is linear in the 
density, we may use the average value for the density as long 
as we make sure that the total path traversed while the velocity 
of the star is small (the contribution to the integral is negligi- 
ble otherwise) is long compared to the average bubble diameter. 
This is indeed the case for an interesting part of the parameter 
space in Vaz,o/Vcirc: We parametrise the accretion suppression by 
/ =< Mflondi > /Meondi and show the path-length per orbit per 
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Differential orbit length over accretion suppression 
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Fig. 4. Differential path length over accretion suppression fac- 
tor / =< Mflondi > /^Bondi- The five curves labelled by their re- 
spective Vaz.o/^circ" value show the path length in parsec that 
a star spends at a given accretion suppression factor / per or- 
bital period and per decade in log(/). The range of /-values that 
our model requires is to the right of the solid vertical bar. The 
horizontal dotted lines show the limiting path length for statis- 
tically uniform accretion. The plot shows that the accumulated 
path length is sufficient for uniform accretion for all orbits with 
I'az.o /I'circ S 3%. For orbits with smaller values of Vaz.o/i'drc ac- 
cretion varies strongly from star to star, and the average value 
shown in Fig.|5]would only apply to the ensemble. The plot also 
shows, that for all orbits with Vaz.o/vcirc < 10% there is always a 
basic amount of gas that is uniformly accreted on to each FRMS 
disc. This amount is sufficient for our model. The part of the ac- 
creted gas that varies from star to star comes only on top of this 
basic amount. See section [J.4.3l for details. 



decade of log(/) in Fig. |4] The orbital period for stars which 
have the outer turning point near the core radius is 0.16 Myr 
for all eccentricities. The minimum accretion time is given by 
the equatorial ejection phase of the 120 star, 2.45 Myr. The 
maximum accretion time would be 4.6 Myr, which would cor- 
respond to the 40 Mq star and the case when the massive stars 
turn into black holes without energetic SNe. Hence, the FRMS 
complete 15 to 29 orbits in the GC during the time when they 
have equatorial ejections. In order to accumulate a total path at 
low velocity of a bubble diameter, they must therefore spend at 
least 0.008 pc or, respectively, 0.015 pc at low velocities, per 
orbit. These two limits are shown as dotted horizontal lines in 
Fig. H A vertical solid line marks the minimum accretion rate 
we must demand to explain the relative contributions of pris- 
tine and ejected gas in the presently observed GC stars. It can 
be seen that, for small values of Vaz.o/i'ciic. the differential path 
length distribution is well described by a power law. 



di 



dlog(/) 



-f 



-3/5 



with an upturn at high values of log(/). Figure |4] shows that 
the accumulated path is generally above the statistical limit for 
Vaz.o/Vcirc ^ 3%. For such orbits, Eq. [T2]is therefore valid for 
each star. Below this value, Eq. [T2]is only valid for an ensem- 
ble that is sufficiently large. Figure|4]also shows that for all stars 
with Vaz,o/i'ciic < 10%, the integrated path between log(/) = -3.5 
and log(/) = -2.5 is always much larger than the statistical path 



Fig. 5. Dependence of the time averaged accretion rate on to the 
FRMS as a function of the ratio of the azimuthal velocity at the 
outer turning point, Vaz.o, over the circular velocity Vcirc, for the 
60Mo star (see section|3]4]for details). The solid line shows the 
elliptical approximation. The pluses show the result of a direct 
orbit integration in the actual Plummer potential. There is excel- 
lent agreement between the two methods. 



length limit. Therefore, such orbits have a uniform basic accre- 
tion rate, which is sufficient in the context of our model. Only on 
top of this, there can be a strongly varying component. Within 
these limitations, we may therefore use Eq. [12] to describe the 
average accretion rate. 

The orbit-averaged accretion rate for the 60 star as a 
function of Vaz.o/i'circ is shown in Fig. |5] We have also derived 
the orbit-averaged accretion rate by direct numerical integration 
of the orbit in the Plummer potential, which is shown as pluses 
in the same plot. The good agreement is expected, as the inner 
regions, where the elliptical approximation is bad, do not con- 
tribute much to the average accretion rate. We note that the time 
averaged accretion rate may be well approximated by: 



Mfiondi = 0.7 



Vaz.0/Vc: 

0.1 



MqI Myr 



M 



0.1/\60Mo/ \106mpcm 



lOOK 



3/2 



(13) 



which is fitted to the curve in Fig.|5] 

While the characteristic velocity vq leads to a suppression of 
the accretion rate by about six orders of magnitude, high orbit 
eccentricity may compensate this to some degree. An accretion 
rate of about 1.5 Mo/Myr, as required to get the right dilution 
(i.e., that requested to explain the observed Li-Na anticorrela- 
tion) of ejected and pristine gas (compare Eq. dS)), is reached 
for Vaz.o / Vciic =7%, for the 60 star. In this case, the gas in the 
disc around the 60 star would consist - on average - of 74 % 
ejected gas from the parent star, and 26 % of accreted pristine 
gas. The accretion rate grows quadratically with the mass of the 
accreting star, whereas for the equatorial mass ejection rate the 
dependency is linear. The pristine gas fraction varies therefore 
with the mass of the parent star. At 60 Mq it is close to minimal, 
for 40 (120) Mq it rises to 35 % (31 %), respectively, assuming 
again Vaz.o/vcirc =7%, We discuss this further in section U.6.21 
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3.5. Viscous merging of accretion and decretion discs 

Because of the angular momentum that is likely present in 
the gas, it will settle into the disc at some radius r, which is 
likely large compared to the typical locus of the equatorial stel- 
lar ejecta. Viscous processes are then responsible for transport- 
ing the material within the accretion disc. For a steady, thin, 
Keple rian disc, the viscous timescale is given by dFrank et alJ 
iOOlpSS): 

VGMr 




where a typical value for o- is 0. 1 and we have taken a radius r 
of order the average distance between massive stars. 

Thus, the viscous evolution of the disc is happening on a 
similar timescale as the equatorial mass ejection. The re-formed 
discs near the massive stars are fed by both, material nuclearly 
processed in the massive stars and expelled via the equatorial 
mass ejections, and by pristine gas via accretion in the shad- 
owed equatorial region. It is thus plausible that such discs de- 
velop a mixture of gas that has overall a very comparable contri- 
bution from both, pristine and processed components, similar to 
the composition of the second gene ration low-mass star s, which 
are presently observed in GCs (e.g. lPrantzos et al1l2007b . 

3.6. Star formation in fhe discs 

Discs that are comparable in mass to the central object are known 
as self-gravitati ng discs (compare e.g. Sect. 3 .1 in the re view by 
lArmitagj|201 lb . The Toomre criterion (e.g. IShulll992l) readily 
shows that the inner discs are expected to be gravitationally un- 
stable: 

^^"°'^(lOkms-')(lOMo) (eOMo) (lOO/?o) ' ^^^^ 

where m is the disc mass, which we assume to be equally dis- 
tributed within the scale radius r for this estimate. We have as- 
sumed a higher sound speed in the disc than calculated for the 
gas on larger scales, because near the massive stars, the radi- 
ation will heat the discs (compare e.g. [Xhmed & Sigut 2012). 
However, even for a sound speed of 10 km/s, the inner discs 
reach the critical mass for gravitational instability on a timescale 
of 10^ years. 

3.6.1. Consequences of gravitational instability 

Gravitationally stable, as well as unstable discs have been 
studied in the context of planet formation. The latter case is 
ho wever not yet fully understood (compare e.g. the review 
by iKlev & Nelson|[2()12) . Spiral density waves with associated 
torques and strong accretion are a likely outcome. It is therefore 
possible that some disc material is accreted on to the FRMS, par- 
ticularly in the "off-phases" of the equatorial ejections. Yet, this 
material should then also have a high residual angular momen- 
tum compared to the gas in the outer regions of the star. Thus the 
time interval to the next ejection should be shortened. However, 
the time interval between the ejections is already much shorter 
than the viscous timescale ( Eq. fT4l i and the timescale for the 
disc to become gravitationally unstable. Hence, the disc sees an 
almost constant flux of material from the star. It is thus possible 



that the main sink for the disc mass would be the formation of 
gravitationally bound objects. 

Likely, the disc will remain near ~ 1, which implies a 
disc mass that declines with time, because also the FRMS looses 
a large fraction of its mass. If the second generation stars would 
form in this way, one would therefore expect that the formation 
rate of second generation stars, or proto-stars which may con- 
tinue to accrete for some time, is highest towards the beginning 
of the equatorial ejections, i.e. about 2-3 Myrs after the birth 
of the first generation of stars (1 Myr till the equatorial ejec- 
tions start also in the lower mass FRMS, a; 1 - 2 Myr viscous 
timescale and timescale for the disc to become gravitational un- 
stable). The equatorial ejections cease entirely about 6 Myr after 
the birth of the first generation of stars, i.e. 2 Myrs into the SNe 
phase, if massive stars explode (compare Sect. |4}. We show in 
Sect. 14.41 below that turbulence in this phase stops accretion on 
to the discs. This effect also keeps the second generation stars 
chemically clean from supernova ejecta. Once the discs are no 
longer fed, they are expected to form their gravitati onally bound 
objects and disappear within at most a few Myrs (iHaisch et al.l 
2001). The formation of the second generation stars is therefore 
complete at, roughly, 10 Myrs after the birth of the first genera- 
tion of stars. 

Thus, while we are not in a position to add all the details 
that still need to be worked out in the context of the formation of 
self gravitating bodies in discs, it appears viable that the second 
generation of stars forms over a timescale of 6-10 Myrs in discs 
around FRMS (Fig.[Tl second row from top). 

3.6.2. The distribution of orbit eccentricities and the ratio of 
pristine-to-processed gas in the second generation 
stars 

We have seen in Sect. 13.4.31 above, that the amount of gas a 
given FRMS disc accretes from the ICM depends strongly on 
the parameter Vaz,o/i'circ ■ The distribution function of Vaz.o/i'circ 
for the sample of FRMS in the cluster will therefore directly 
determine the distribution of the amounts of accreted ICM to the 
FRMS discs. 

The actual abundance pattern in each second generation star 
will however not only be determined by the overall ratio of pris- 
tine gas accreted from the ICM over ejected gas from the FRMS. 
The abundances in each disc will be also a function of time be- 
cause of the time dependence o f the abundances in the equatorial 
ejecta dPecressin et al.ll2007bl) . The composition of the equato- 
rial ejecta is first very similar to the one of the pristine gas. As 
time goes on, the composition changes and shows incr easing 
signatures of hydrogen burning. iDecressin et al. IdlOOTb) have 
shown that the observed spread in light element abundances 
may only be obtained if the time variation of the abundances 
in the ejecta is used. For example, only at late times do we 
yield oxygen depletions high enough to be compatible with the 
oxygen poor end of the abundance distributions in many GCs 
(lGrattonetal.1120121) . If we would mix all the ejecta of a given 
star into a common reservoir, which would be the disc in our 
case, the average oxygen depletion would already be insufficient 
to explain many of the observed second generation stars. Further 
mixing with pristine gas would make the situation even worse. 
We therefore require that the star formation in the FRMS discs 
happens sequentially. This means that for each second genera- 
tion star, its formation is complete on a timescale significantly 
shorter than the lifetime of its parent FRMS disc, so that it may 
conserve the current abundance pattern of the disc at the time 
of its formation. This may for example be achieved via migra- 
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tion towards outer orbits, where the gas densities in the discs are 
low, or encounters between second generation stars which might 
scatter them out of the plane of the disc. 

A consequence of this way of star formation is that the dis- 
tribution of abundances will always be continuous. Even in the 
extreme case when we assume a population of stars with cir- 
cular orbits, which therefore do not accrete anything from the 
ICM, together with a population of stars with low eccentrici- 
ties, which therefore accrete strongly, we would only expect that 
for the population which have more ICM accretion, the range of 
metallicities would be reduced. It is difficult to imagine that this 
co uld lead to a double p eaked distribution like the one observed 
by iMarino et al.l (l2008h in M4. We could however well explain 
the m ore uniform distribution found in NGC 1851 (iLardo et al.l 
I2OI2I) . Details of abundance distributions will be the subject of 
future work. 

4. Supernovae 

It is well established th at massive s tars (M > 25Mq) end their 
lives as black holes (e.g. iFrye ri 19991). Stars initially less massive 
than this are expected to produce energetic S Ne-events or, may 
be. Gamma ray bu rsts for fast rotators (e.g. lYoon et all l2006t 
iDessart et al.ll2012}). This is much less clear for t he more mas- 
sive stars (e.g. lFrye rl ll999llBelczvnski et al.ll20T2l) . which might 
turn silently into bla ck holes (compare also the discussion in 
iDecressin et al.ll20Tol and references therein). 
This introduces a considerable uncertainty into our model time- 
line: The 120 Mq star would explode at 3.5 Myr, the 25 Mq 
star at 8.8 Myr. The wind-bubble phase ends and the supernova 
phase begins at some point within this range with the first ener- 
getic supernova. In the following we will use the term supernova 
for simplicity, always subsuming the possibility of Gamma-ray 
bursts. 

We will now first argue that the SNe will likely not signifi- 
cantly alter the picture of the star formation in the FRMS discs 
(Sect. 14. iV The SNe will cause substantial turbulence in the ICM 
(Sect. 14. 21 ). which will likely lead to mixing of the SN-ejecta with 
the cold phase of the ICM (Sect.l43Tl. The SN-ejecta may how- 
ever not enter the second generation stars, because on the one 
hand, the ICM does not form stars in this phase, and on the other 
hand, the ICM can no longer accrete on to the FRMS discs be- 
cause of the large gas velocities (Sect. l4!4] i. From the end stages 
in the lifes of the FRMS onwards, we expect the second gen- 
eration stars to be dispersed, first within the half-mass radius, 
and on the relaxation timescale also throughout the whole clus- 
ter (Sect. |4j~ 



SNe driven turbulence 



4.1. The effect of supernovae on the associated star-forming 
disc 

Supernova ejecta are fast and cannot be retained in the grav- 
itati onal potential of the di scs around their parent stars (com- 
pare IDecressin et al. I l20073bh . Since the disc occupies only a 
small solid angle, much of the energy will escape into the 
other directions, especially if the energy release occurs in the 
form of iets as might be the case for fast rotating massive 
stars ( e.g. iBisnovatv -Kogan & Moiseenkol2()0^lTakiwaki et all 



3gai 

l2009t iDessart et al. 20121) . and ablate onlv the surface of the 
disc. At this point (4.5 Myrs for the 6OM0 star) we might expect 
that a large fraction of the inner disc gas has already accreted on 
to the newly formed stars. In this case, some remaining debris 
might be cleared away by the explosion. Yet, it is also possible 
that a substantial amount of gas is left. If there would still be a 





DU 


CO 


40 


E 








> 


20 









0.9 




Q. 


0.6 








0.3 




0.0 




10 20 
Time / Myr 



30 



Fig. 6. RMS-velocity (upper part) and characteristic radius 
(lower part) established in the SNe driven convective phase. The 
left axis shows, respectively, km/s and pc. The right axis indi- 
cates, respectively, velocity in units of the central escape velocity 
and radius in units of the half-mass radius. The solid black line is 
for the case when all massive stars produce energetic SNe. The 
dotted blue line assumes direct black hole formation without en- 
ergy release to the ICM for all stars with M > 25Mq. The plot 
shows that in the latter case, the onset of turbulence is delayed 
by 5.3 Myrs. See Sect. l4.2l for more details. 



substantial amount of g as in the disc, one w ould expect, similar 
as in the simulations o f lGaibleretai](l2012l) for galactic scales, 
a compression of the disc and stronger gravitational instability 
during this phase. 

4.2. The effect of supernovae on the general ICM: turbulence 

Supernovae are energetic enough to overcome the gravitational 
potential of the cluster. They are exp ected to f orm a superbubble 
with an expanding shell. As shown in lPaper A such a shell would 
be Rayleigh-Taylor unstable. The situation is simi lar to the one 
describ ed by the 2D-axisymmetric simulations of iKritsuk et al.l 
(120011) : The fragments of the shell are expected to fall back into 
the GC, and a convective cooling core is expected to form. Now, 
the hot high entropy gas escapes, while cooling shell fragments 
are continuously formed at the edges of escaping bubbles, and 
fall back towards the centre, much like in a boiling pot of water. 

We will now estimate, if the typical turbulent velocities could 
lead to a significant increase of the scale height of the gas. 
We are particularly interested in the cold gas, since most of 
the mass is expected to be in this phase. Hydrodynamic en- 
ergy transfer is expected between the different gas phases (e.g. 
Krause & Alexander 20(33)- In equilibrium, the energy input by 
SNe should be balanced by radiative losses. But it is not clear 
th at such an equilibriu m may be reached: In the simulations 
of iKritsuk et al.l (l200lh the convection z one expands slowly, re- 
flecting a steady energy accumulation. iKrumhoiz et all (120061) 
summarise simulations of isothermal turbulence, and find a de- 
cay time of fdis - 0.83/l,„/v'rms- The injection scale in our case 
is given by the half-mass radius, since this is roughly the scale 
where the she lls are destroyed by th e Rayleigh-Taylor instabil- 
ity. Following lKrumholz et al.l (120061) . we set the injection wave- 
length to Ain - 4ri/2, because the gas is first in isotropic infall, 
and also because this is the maximum possible wavelength. The 
rms-velocity is given by £ = (1 - esf)MtotV^s/2, because the cold 
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gas is supersonic so that the thermal energy may be neglected. 
We may thus approximate the change of the energy in the gas by 



dt 



(16) 



Once the rms-velocity is calculated by this approach, the char- 
acteristic radius out to which a given gas filament may ascend in 
the potential of the GC follows from: 



^ rms 



2 

e,0 



(17) 



where the central escape velocity Ve,o = 183 km/s for our model 
GC. We plot the evolution of the rms-velocity and the charac- 
teristic radius of the convective flow, for the case when the gas 
would be supported by turbulence alone, in Fig.|6] Our choice of 
Ain guarantees the smallest possible dissipation rate. The turbu- 
lent velocities reached are large compared to normal interstellar 
medium (ISM) turbulence, about 50 km/s. Yet, this falls short 
of the escape speed. However, if the gas would have some other 
means of support, apart from ram pressure due to turbulence, 
it would be quite likely, that it would still be supported in this 
other way. Alternative means of support against gravity could 
be via magnetic fields or to some degree also radiation pres- 
sure dKrumholz & Thompson 2012). If this would be the case, 
our result might also be interpreted in the way that SNe driven 
turbulence is not able to increase the scale height of the gas in 
this phase in any significant way. We note that the characteris- 
tic turbulent velocity is mainly determined by the ratio of en- 
ergy injection to gas mass, and also by the size of the cluster via 
fdis- Thus, since GCs all have a typical radius of a few pc, and 
since the ratio of energy injection to gas mass depends only on 
the star formation efficiency, which should be al ways around 1/3 
for effi cient ejection of the first generation stars (iDecressin et al.l 
1201 Ol) . we expect essentially always turbulence at 50 km/s. For 
proto-cluster clouds with initial masses above IO^Mq, which 
should be approximately the minimum initial gas mass to form 
a GC, this is below the escape velocity. Thus, apart from very 
low-mass protocluster clouds with high star formation efficiency, 
these findings should be generally applicable for all GCs. 

4.3. Mixing of supernovae ejecta witli pristine gas 

Turbulence is consid ered to be a decisive factor in mixing theo- 
ries (for a review see IScalo & Elmegreenll2004l) . Mixing of stel- 
lar ejecta into ICM that later forms subsequent generations of 
stars is strongly constrained by observations (e.g jGratton et al.l 
1201 2l and sect. [1] above). Here, we are interested in an order 
of magnitude estim ate and therefore use the mixing length ap- 
proach. We follow IXie et al.] ( Il995h who apply mixing length 
theory to interstellar clouds. Within this framework, the mixing 
timescale is given hy t x H/ Vd, where H represents a mean of 
abundance and density scale heights, and is the diff'usion ve- 
locity, given by Vd ~ K/H. The diffusion coefficient K is simply 
approximated by the product of characteristic turbulent veloc- 
ity, Vt and the mixing length L. Let us assume filament thick- 
nesses of order 0. 1 pc, in the SNe driven turbulence phase, and 
set L = // = 0. 1 pc. Thus, the diffusion velocity equals the turbu- 
lent one. We have shown above that the expected turbulent veloc- 
ity in the SNe driven turbulence phase is approximately 50 km/s, 
which therefore would be the characteristic diffusion velocity in 
the cold gas. The sound speed in the hot component wi ll be larger 
than this. X-ray observations (e.g. Ijaskot et"ani201 ll) generally 



find hot gas temperatures of order 10^ K in superbubbles. The 
sound speed will consequently be several 100 km/s. Cold gas 
will therefore more efficiently be mixed into the hot gas than vice 
versa. Using these assumptions, the mixing timescale is given 
by: 



:2000yrs 



H 



0.1pc/\50km/s 



(18) 



This should be compared to the dynamical timescale, i.e. the 
time needed by bubbles enriched with fast winds and SNe ejecta 
to escape the GC. At times when the GC is highly turbulent the 
dynamical timescale will be given by: 



T"dyn,t = n/l/Vt 

~ 60,000yrs 



In/2 



\3pc / \50km/s 



(19) 



When the turbulent velocities are small, the gravitational accel- 
eration on the cold gas which needs to refill the bubble vol- 
ume will be the limiting factor Therefore, the Rayleigh-Taylor 
timescale sets the minimum for the dynamical timescale: 



Tdyn,RT = ^j^^njlfg 



60,000yrs(^ 



1/2 



\3pc/ \5x 10-*'cm/s' 



-1/2 



(20) 



where we have scaled the gravitational acceleration g to the av- 
erage value within the half-mass radius for our model cluster 
Thus, the dynamical timescale is always around 60,000 yrs. 

This means that in the supernova phase, we expect the SNe 
ejecta to mix with the cold gas phase of the ICM, because of 
fast turbulent mixing. On the other hand, in the absence of an 
efficient driver, any turbulence should decay quickly and the gas 
velocities are expected to be of order the sound speed, which is 
expected to be about 1 km/s (compare Sect. I3.4l i. The mixing 
timescale increases therefore to about 100,000 yrs, almost twice 
the dynamical timescale. Hence, we expect no significant mixing 
in the wind phase before the onset of the SNe. 

We note that the issue of turbulent mixing in the ISMcan- 
not be regarded as settled yet. De Avillez & Mac Low (l2002h 
find in hydrodynamic simulations of SNe-driven turbulence in 
the Milky Way disk diffusion coefficients which are one to two 
orders of magnitude smaller than in the mixing length approach. 
Correspondingly, the mixing timescales we give above would 
underestimate the true timescales by one or even two orders of 
magnitude. While an increase of one order of magnitude would 
not change our result, two orders of magnitude would mean that 
also in the SNe-driven turbulence phase, no significant mixing 
would occur 

It has been suggested that stars with masses above 25Mo 
might not form SNe, but would instead collapse to a black hole 
without noteworthy energy release (compare above). This would 
delay the onset of strong turbulence until about 8.8 Myr after 
the birth of the first generation stars, the time when the 25 
stars would explode. The FRMS could in this case accrete from 
the ICM up to this point. Yet, the lack of turbulence would also 
prevent mixing, so that the stellar ejecta carried by the fast ra- 
diative winds of the more massive stars at this time, which carry 
He-burning products that are not found in the second generation 
stars, cannot enter the FRMS-discs. The same argument applies 
for the end of the wind bubble phase for the case when the high 
mass FRMS do develop SNe. 
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4.4. The effect of turbulence on accretion 

As discussed in Sect. 13.41 accretion is significantly suppressed, 
if the relative velocities between the accreting objects and the 
accreted gas are large. This is the case for the SNe phase, as 
already the typical gas velocities are around 50 km/s. Thus, ac- 
cretion on to the FRMS discs may only happen during the first 
few Myr (4 to 9 Myr, depending on the mass limit for energetic 
SNe events) after birth of the first generation of stars, before 
the supernova phase. The associated smaller amount of accreted 
gas may however be compensated by a slight increase of the or- 
bital eccentricities of the FRMS (Sect. 13.41 ). We have shown in 
Sect. l3.6l that the FRMS discs have disappeared long before the 
end of the SNe phase. Once a given FRMS has exploded as a 
supernova, there is therefore at most a short period of activity, 
and the dark remnants then remain in-active until the last SN has 
exploded. 



4.5. Dispersal of the second generation stars 

The FRMS loose about half of their mass via equatorial ejec- 
tions loaded with hydrogen-burning products during the main 
sequence and the luminous blue variable stage, and this is the 
main mass loss channel in this early phase. We expect that to- 
wards its end, each FRMS disc has about the same mass as the 
FRMS, perhaps even a bit more. A considerable fraction of the 
disc mass should already have been used up to form the sec- 
ond generation stars. Thus, even towards the end of the equato- 
rial ejections, it is expected that gravitational interactions lead to 
migration of the second generation stars in the discs, and quite 
possibly also to some ejections from the discs. 

From the end of the equatorial ejection phase, the FRMS go 
on loosing mass via fast radiative winds and perhaps eventual 
SNe ejections. The final dark remnants might then have only 
about ten per cent of the mass that the FRMS had at the end 
of the equatorial ejection phase. Consequently, the gravitational 
binding energy, which is to first order simply proportional to the 
square of the total mass, would have decreased by a factor of a 
few. Virial equilibrium at the end of the equatorial ejection phase 
demands that the kinetic energy is half of the binding energy. 
Therefore, after a SN-ejection, the total energy of the system 
would be positive. Even without SN-ejection, the mass loss of 
the FRMS after the equatorial ejections have ceased is consider- 
able. Consequently, the second generation star orbits widen. This 
increases the cross sections for collisions. Most of the second 
generation stars may be lost from the parent star, but would re- 
main centrally concentrated within the half-mass radius, as also 
the FRMS and the dark remnants (compare Fig. [TJ. The mixing 
of the first and second generation stars due to two-body relax - 
ation eff'ects has been studied in detail bv lDecressin et al.l(l2008l) . 
They show that the two stellar populations need two relaxation 
timescales to mix, which corresponds to about 2.8 Gyr, from 
their Eq. (1). Therefore, there is not enough time for the second 
generation stars to spread over the cluster, before the majority 
of the first generation stars are lost from the cluster during the 
gas expulsion phase (section |5]l. Yet it appears also well possi- 
ble that the dark remnants retain some of the second generation 
stars, which might thus form X-ray binaries later. X-ray binaries 
are well known in GCs (e. g. iBamard etal1l2012t lBarretll2012t 
IStacev et al.ll2()T2h . 



5. Accretion on the dark remnants 

After the last type II SNe has exploded, roughly 35 Myr after the 
formation of the first generation of stars, turbulence decays. With 
the formula given in Sect. l4.2l above. the decay time evaluates to 
0.2 Myr/v5o, where V50 is the rms-velocity in units of 50 km/s. 

In principle, black holes and ne utron stars may rece ive na- 
tal kicks as large as ^ 400 km s~' (iRepetto et al.ll2012h . large 
enough for them to escape from any given GC. H owever, many 
neutron stars are actually observed in GCs (e. g. iBamard et al.l 
I2OI2I: lBarretll2012l: IStacev et al.ll2012h . The recen t detection of 
two b lack holes in the Milky Way GC M 22 (IStrader et al.l 
l2012h has called into question if all black holes receive large 
natal kicks. Especially if black holes form directly without SN- 
ejection, no significant natal kicks are expected. 

In the following, we assume therefore that the dark remnants 
have the same orbits as the FRMS have had before. Therefore, 
accretion should again set in. Assuming 1 .4 Mq for the neutron 
stars, and 3 M© for the stellar mass black holes and accretion 
over the full solid angle, Eq.|9]predicts about 10-100 Mo/Myr. 
If we assume a similar suppression as in Sect. 13.41 above, we 
end up with 0.. 003-0.03 Mo/Myr. The dark remnants will start 
their energy output after about a viscous timescale, which is 
needed by the gas to get close enough to their surface, or re- 
spectively, the horizon. The viscous timescale is about 1 Myr, 
according to Eq. [14] Therefore, each dark remnant will receive 
0.. 003-0. 03 Mq of gas and then get activated. This compares to 
the Eddington accretion rate of 0.02MDRMo/Myr, where Mqr 
is the mass of the dark remnant in solar masses. This means 
that once loaded with gas, the dark remnants could be active 
at the Eddington rate for about 0.1-0.5 Myrs (Fig. [1] bottom). 
We have shown in lPaper j that such a sudden onset of accretion 
on the stellar mass black holes and neutron stars may plausi- 
bly release a sufficient amount of energy to eject the ICM in at 
most 0.06 Myrs, corresponding to about half a crossing time. 
The mass accreted on to the dark remnants may thus plausibly 
be sufficient for this to occur. The expulsion timescale is short 
enough to remove also the outer first generation low-mass stars, 
which are not to appear in the GCs today. We have also shown 
that a gradual activation of the dark remnants as they become 
available, reinforcing the SNe in that phase, is still insufficient 
to overcome the Rayleigh-Taylor insta bility. Fo r the efficiency 
of energy transfer we have adopted in iPaper II (20 per cent of 
the Eddington luminosity), the limiting proto-cluster cloud mass 
up to which the gas expulsion worked was about 2 x 10^ Mq. 
This would mean that gas expulsion would work for most of 
the observed clusters, apart from the most massive ones, e.g. w 
Centauri and M 22, where al so a spread in iron abundances is 
observed dGratton et al.ll2012l and references therein). The re- 
peated star-bursts expected in this case are similar to the predic- 
tio ns for nuclear star clusters without supermassive black holes 
bv lLoose et al. l (fT98l . 

For GCs below this mass limit, further gas accretion events 
are expected from about 40 Myrs after th e birth of ffie first gen - 
eration stars due to the slow AGB -winds (iD'Ercole et al.ll2008h . 
Any such events would re-activate the dark-remnants. Both, the 
accretion rate and the binding energy are approximately linear in 
the gas density , if the stellar potential dominates. Otherwise, the 
binding energy drops even stronger with gas density. Therefore, 
we expect again quick gas expulsion. This would be similar to 
the case of super-massive black holes in nuclear star clusters and 
ellipt ical galaxies, which may also be activated by AGB-winds 
(e.g. iGaibler etal.]l2005t iDavies et all 120071: ISchartmann et al.1 
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l2010t iDavies et alJ l2010l). Black-h ole feedback may suppress 
star formation (e.g. lDi Matteo et al.ll2005h . 



6. Discussion 

6.1. Orbit eccentricities of tlie FRMS 

The weakest point of the presented scenario is perhaps the fine- 
tuning required for the orbital eccentricities of the FRMS, in or- 
der to obtain high enough accretion rates, both, on to the FRMS 
discs and also the dark remnants (Sect. 13.4) . We require veloc- 
ities near the outer turning points of less then about a tenth of 
the circular velocity there. The former evaluate to a few km/s. 
An important point to make in this context is that the require- 
ments are the very same as for the dark remnant accretion, as the 
FRMS turn into the dark remnants at the end of their lives. If gas 
expulsio n is inde ed due to the dark remnants, as we have pro- 
posed in iPaper iL then the orbits are also ok for accretion on to 
the FRMS, unless they were significantly affected by the trans- 
formation (compare sect. |5]above). 

6.2. Support of tiie gas and relation to GC formation 
scenarios 

The FRMS orbits will of course be related to the formation 
of the GC. This is so far u nclear, and different scenarios are 



considered (compare e.g. [Harris et al 



Kravtsov & GnedinI 2005 : GnedinI 2011 : Gray & Scannapiecol 



[T99i iKrausel I2OO2I 



[Harris & HarrisI jj). In general, one may consider two 
classes of physical conditions: First, the first generation stars 
might be formed while the gas is in free collapse. This is 
among the s cenarios currently discussed for star formation in 
general (e.g. IZamora-Aviles et al] l2012h . Then, by definition, 
the stars will be bom with large radial velocities, and much 
smaller azimuthal velocities, roughly as required. The second 
option is to have internal support against gravity in the proto- 
cluster cloud. As the decay time for turbulence is short (compare 
Sect. I4.2I 1 and an unusually high level of turbulence is required 
(x! 100 km/s), it seems unlikely that such a protocluster would 
be supported by turbulence. If it was, the stars would be ex- 
pected to form with comparable radial and azimuthal velocities, 
which would be a serious problem for this scenario. Thermal 
pressure seems also unlikely, as the cooling times would be 
very short. The protocluster cloud might also be supported by 
magnetic fields. In general ISM simulations, the coldest and 
hence densest phase is found to be magnetically supported, 
in contrast to turbulent support for intermediate temperatures 
dde Avillez & Breitschw erdt 20 (j5|). This a grees with observa- 
tions of molecular clouds (lCrutcheijl2012l) . In a magnetically 
supported protocluster cloud, the stars would lose their magnetic 
support immediately after formation, and one would also expect 
highly eccentric orbits. 

6.3. Effects of magnetic fields 

Magnetic support would mean that the proto-cluster cloud could 
keep the same scale height as the stars until the gas is expelled 
from the GC. This would be highly beneficial, as for our scenario 
accretion takes place mainly when the stars are in the vicinity 
of the core radius. Magnetic fields might however also inhibit 
accretion. iCun ningham et al. (2012) show that Bondi accretion 
from an isothermal, magnetised gas is suppressed by a factor 
0.4cs/ca, where ca refers to the Alfven speed. This might seem 



moderate compared to the (cs/v)^-factor for bulk motions at ve- 
locity V. Yet, to support the protocluster cloud magnetically, this 
also amounts to two orders of magnitudes, which would signif- 
icantly restrict the available orbit eccentricities. Ambipolar dif- 
fusion should however be important in this context, as the ma- 
j ority of the cold gas would be neutral. Also. ICunningham et al.l 
(1201 2h assume that the magnetic flux accumulates around the 
accretor, which is the underlying physical process inhibiting 
the accretion. In our case, we would expect the flux to be ad- 
vected into the accretion discs, there to be dissipated by tur- 
bulence. Qualitatively, we therefore believe that a magnetically 
supported protocluster cloud would be a viable alternative in 
the present context. Radiation pressure might also help to sup- 
port the gas after the first gene ration massive stars have formed 
(iKrumholz & Thompsonll2012h . 



6.4. Stronger mass segregation as alternative to supporting 
the gas on a scale of the half-mass radius 

An alternative to supporting the gas on the scale of the half- 
mass radius would be to assume a stronger concentration of 
the gas as well as the massive stars towards the centre of the 
GC, implying stronger mass segregation: In our model setup 
(Sect. |5]), we specified that the massive stars are confined to the 
sphere inside the half-mass radius. In fact, iLeigh et aP (l20 1 3h 
in their model for dark-remnant accretion and dynamical black- 
hole ejection in GCs, assume a much stronger mass segregation. 
This might happen, if a GC is formed from merging sub-clumps. 
In this case, strong m ass segregation of the stars is expected to 
arise in about lO *^ yr (lAlUson et al.ll2010l:lGirichidis et al.ll2012l: 
iPang etalJl2013h . Assuming a velocity dispersion inversely pro- 
portional to the stellar mass, as in Leigh et al. (2013), would re- 
sult in our relevant FRMS population having maximum veloci- 
ties of order the sound speed. In such a scenario, both, the mas- 
sive stars and the gas could be very much concentrated towards 
the cluster centre. Hence, accretion would no longer be affected 
by the orbital parameters, as the FRMS velocities would always 
be small. Apart from the stronger accretion, the formation of the 
second-generation stars would in this variant of the scenario pro- 
ceed in almost the same way as in the standard case (Sects. 13.51 
and l3.6l l. The star forming discs must have some way to prevent 
complete mixing of the ejected and the accreted gas in order to 
explain the most oxygen poor stars (compare Sect. 13. 6. 2) . The 
SN-phase would be essentially unaffected, whereas the accretion 
on to the dark remnants would also be enhanced. The latter might 
lead to an observationally interesting collective active phase of 
order the viscous timescale, i.e. about 1 Myr. 



6.5. Limitations of dark-remnant accretion due to local 
radiative feedback 

Radiative feedback has been found to limit accretion on to 
stellar-mass black holes in spherical and clumpy accretion flows 
(iMilosavlievic et aT]l2009h . This applies to the dark-remnant ac- 
cretion phase, as we have assumed spherical accretion in Sect.|5] 
Milosavlievic et al. (2009) show that the average accretion rate 
may be limited by a value two or three orders of magnitude 
below the Eddington limit. They caution however that angular 
momentum, turbulence, and thermodynamic phase segregation 
might affect the result. They also find that on short timescales 
the accretion rate may exceed their limit substantially. We re- 
quire only a very small active phase of the dark remnants to ex- 
pel the gas (< 10^ yr). Further work is clearly required to settle 
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the case. The strong observational constraints on the present- 
day gas content of GCs and the Fe-uniformity of the present- 
day stars put stron g constraints on the gas expulsion mechanism 



aay stars nut stron g 
(compare iPaper J) . 



6.6. Top-heavy IMF 

Gas expulsion could be substantially eased by a non-standard, 
top-heavy IMF. For the model presented in this paper (Sect. [Sj, 
we have used a normal IMF, following oth er recent work (e.g. 
iDecressin etal] 120101: iD'Ercole et "all 1201 ll) . However, the first 
generation low-mass stars are of no importance for our scenario. 
In fact, there is even an issue with their ejection, which we have 
found here not to work with SN-feedback alone (compare pa- 
per I and Sect. |5] above), and which may therefore require the 
coherent activity of the dark remnants. The ejected low-mass 
stars might contri bute significantly to the stella r population of 
the Galactic halo ( ISchaerer & Charbonnell201l1) . but there is no 
strong constraint on how many stars have to escape. The IMF for 
the first generation stars might the refore have been top-heavy 
(e.g. iPrantzos & Charbonnell 120061 and references in Sect. [1}. 
This would not necessarily be in contradiction with recent work 
on cl usters formed under typical present-day M ilky Way condi- 
tions ( lHennebellel2012l:lKrumholz et alj2012l) . as the conditions 
in which the Galactic -halo GCs formed are not known. The cen- 
tral disc in the Milky Way's nuclear star cluster, a massive and 
concentrated star cluster similar in this respect to GCs, is an ex- 
ample where observ ational evidence for a top-heavy IMF exists 
dBartko et al.ll201(3h . A high star-formation efficienc y, as often 
assumed for GCs, is also discussed for this region (e.g. lSiUc et al.l 
l2012h . A top-heavy IMF has been suggested in the context of 
GCs on the basis of the energy requirements for the gas expul- 
sion dMarks et al.ll2012l) . 

In the extreme case, assuming that almost no first generation low 
mass stars form, the initial mass of our model-protocluster cloud 
(9 X IO^Mq, Sect. |2]i could be reduced by one order of magni- 
tude. Because of the combined effect of lower gas densities and 
smaller stellar velocities, the accretion suppression factor would 
be constrained to be larger by a factor of a few, moving the verti- 
cal line in Fig.|4]to the right. This would nan^ow the allowed or- 
bital parameters, requiring somewhat higher eccentricities (but 
see Sect 16.41 for a non-exclusive model variant which has the 
opposite effect). Depending on the star-formation efficiency and 
the energy-injection efficiency in the SN-phase, the SNe might 
already eject the gas. We stress that this refers to the extreme 
scenario. We have shown in paper I that even tripling the mas- 
sive stars would not yet suffice to expel the gas via supernova 
feedback. The argument is still valid that favourable conditions 
for accretion on to the FRMS imply later accretion on to the 
dark remnants. A possibly reduced gas density due to successful 
SNe-feedback would cause less accretion on to the dark rem- 
nants (linear dependency, compare Eq. (|9)), but the binding en- 
ergy would also be reduced. Hence, the dark-remnant feedback 
would be expected to keep the GC gas-free. 



6. 7. Comparison to Be-star decretion discs and "first-star" 
simulations 

A key point in our scenario is the formation of the second gen- 
eration stars in discs, fed by both, accretion of pristine gas and 
equatorial ejections of the FRM S. Such equatorial ej ections are 
known from Be-type stars (e.g. iRivinius et al. I I2OOII) . Recently, 
the spreading out of such discs is also attempted to be under- 



stood in terms viscous evolution (iHaubois et al.ll20l"2l) . similarly 
as we propose here for the FRMS discs. 

The formation of stars as opposed to planets in circumstellar 
discs might seem unusual. Yet, detailed models have recently 
be comput ed in the context of the formation of the first massive 
stars (e.g. iGreif et aL 2012h . which also m ay have been very 
fast rotators (IStacv et al.ll20lTl) . For example. IStacv et al.l (120121) 
simulate the formation of a a; 3QMq star with radiative feedback. 
They find that the radiative feedback does not stop the secondary 
star formation in the accretion disc. During their simulation time 
of a few 1000 yrs, th ey form a few pro to-stars of sometimes 
several solar masses. IClark et all (1201 ll) form four secondary 
proto-stars within 110 yrs in a similar simulation. Some of these 
proto-stars may end up in the primary star. Adapted to the FRMS 
case, this might then lead to a more rapid sequence of ejections. 
Clearly, a lot of details need to be worked out, some are not 
comparable (e.g. the first stars are usually assumed to form in 
dark matter halos). Yet, we believe that these developments are 
encouraging. 



7. Conclusions 

Using FRMS data, we have developed a comprehensive sce- 
nario for the formation of second generation GC stars (compare 
Fig.©: 

- The ICM obtains a spongy structure in the wind phase 
(Fig. [1] second row) due to the space-filling wind bubbles 
of the massive stars. 

- We confirm the analysis of Conroy and Spergel that the 
Lyman- Werner photon flux density is sufficiently high that 
most of the hydrogen molecules dissociate and no stars may 
form in the normal way, during the wind and the supernova 
phase. 

- If the first generation massive stars have however equatorial 
ejections as expected in the case of fast rotation, we show 
that accretion on to their discs resumes in the shadow of 
the equatorial ejecta. Massive discs are formed fed by both, 
FRMS ejecta from the inside, and accretion of pristine gas 
from the outside. Within an order of magnitude estimate, 
both contributions can b e comparable, as requir ed by self- 
enrichment calculations (IDecressin et al.ll2007bl) . if the or- 
bits of the FRMS are sufficiently eccentric. 

- The second generation stars may then form due to gravita- 
tional instability in these discs. (Fig. [1] second row). This 
second-generation star formation might carry on for a few 
Myr through the supernova phase, because SNe ejecta are 
not expected to enter the discs. The increased pressure due to 
the SN may however also compress its associated disc such 
that the last stars in the respective disc form at this occasion. 

- The second generation stars form sequentially in their par- 
ent FRMS discs, with the current light element abundances 
of the respective discs. We propose this as the physical 
m echanism underlying the detailed abundance calculations 

ecressin et allboOTbl) . 



of|D 

- The formation of the second generation stars is completed 
latest about 10 Myr after the formation of the first genera- 
tion. The second generation stars are first dispersed within 
the half-mass radius, and much later, after the gas expulsion, 
also throughout the entire GC. 

- The SNe drive turbulence at about 50 km/s (Fig. [Tj third 
row). The remaining gas outside the discs might then mix 
with the SN ejecta due to turbulence caused by the SN en- 
ergy. The gas remains bound to the cluster until the SNe have 
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ceased, turbulence has decayed and the gas can once more 
accrete suddenly on to the dark remnants. The gas may not 
accrete on to the FRMS discs before due to the high turbu- 
lent velocities. In this way, the second generation stars are 
prevented from pollution by SN ejecta and He-burning prod- 
ucts in general. These findings should hold for all globular 
clusters, apart from perhaps the ones at the low mass end if 
their star formation efficiency significantly exceeds 1/3. 

- In the dark remnant accretion phase (Fig. [T] bottom), the 
gas is efficiently expelled due to the strong energy release 
associated with the accretion on to the dark remnants. This 
happens fast enough so that a large fraction of less tightly 
bound first generation stars are also lost. 

Potential problems for the scenario include: 

- For efficient accretion, the orbits of the massive stars need 
to be eccentric, with low angular momenta, and velocities 
near the outer turning points below about 10 % of the cir- 
cular velocity. This may relate to the formation scenario. It 
would pose a serious difficulty for this model, if such low 
angular momenta could be ruled out observationally and, 
consequently, stronger mass segregation would have to be 
assumed (Sect. 16741 1. 

- Further, the gas scale height has to be comparable to the core 
radius, as the main gas accretion would occur near the core 
radius. The necessary suppor t of the gas against gravity may 
be due to radiation pressure (iKrumholz & Thompsonll2012h 
and magnetic fields, details are however beyond the scope of 
this work. Stronger mass segregation would again alleviate 
this problem. 

- Star formation in discs around massive stars is only recently 
being explored. While the results from the "first stars" sim- 
ulations are encouraging, simulations tuned to the specific 
conditions in FRMS are not yet available. 

- If for some reason gas expulsion by dark remnant accretion 
would not work (Sect. 16.5b . all gas should at some point 
form stars from the quite possibly SN enriched gas (Sect l4.3l 
for mixing uncertainties), unless some other way of gas ex- 
pulsion would be found. This would however conflict with 
observations except for the rare and most massive GCs like 
u) Centauri and M22 that exhibit Fe spread. 
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Fig. 1. Sketch of the proposed model of first 40 Myrs of evolution of a two-populations globular cluster together with the global 
timeline (top to bottom) and the associated time when the stellar life ends for selected stars on the right. The leftmost column 
shows the whole cluster. The second column from the left represents a zoom on to a FRMS. Four important stages are depicted. 
First row from top: First a mass segregated star cluster is formed with all the FRMS inside the half-mass radius (dashed line) and 
with the initial gas (light pistachio green shade) remaining after the star formation. Each massive star (the blue interior signifies 
the convective hydrogen burning core, the bright forest green envelope depicts the convective envelope, which has still the pristine 
composition at this stage) creates a hot bubble around it. The corresponding wind shell is shown in a darker shade of green than 
the uncompressed gas. Second row: All the hot bubbles connect and create a spongy-structure in the centre. At the same time, slow 
mechanical winds around the FRMS create a disk around them. The outer envelope of the FRMS is now blue to denote that it has 
been contaminated by hydrogen burning products. In the interaction between the ejected disc (blue) and the accreting interstellar 
medium (ISM, green), a second generation of chemically different stars is born (blue filled circles). Third row: SNe (red stars with 
straight lines) fail to eject the gas but create a highly turbulent convection zone. Further accretion on to the remaining FRMS discs is 
inhibited. The equatorial ejections also end around this time and the formation of the second generation stars is completed by about 
10 Myrs on the global time axis. Bottom row: Later, rapid gas expulsion takes place and removes all the remaining gas together with 
the majority of the less tightly bound first generation stars out of the cluster potential well. Such rapid gas expulsion is likely to be 
due to the activation of black holes by accretion of matter. 15 



